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ABSTRACT 


A  procedure  is  presented  for  the  estimation  of  extreme 
values  of  stationary  Gaussian  random  processes  with 
arbitrary  bandwidths.  This  approach  is  based  on  the 
analytic  envelope  defined  by  the  Hilbert  Transform;  this 
envelope  is  Rayleigh  distributed  regardless  of  bandwidth. 
For  experimentally  derived  data  that  has  been  converted 
into  digital  form,  the  Hilbert  Transform  is  approximated 
using  algorithms  implemented  on  a  digital  computer  to 
produce  samples  of  the  envelope's  time  history.  Next,  the 
degree  of  correlation  between  these  envelope  samples  is 
taken  into  account  using  a  method  developed  from  simulation 
studies  of  a  series  of  synthetic  Gaussian  time  histories 
with  varying  bandwidths.  Once  this  correlation  effect  has 
been  estimated,  the  standard  methods  of  order  statistics 
are  applied  to  these  samples  using  the  Rayleigh  probability 
density  function.  Examples  of  applying  this  procedure  to 
experimentally  derived  data  are  presented. 
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This  work  was  sponsored  by  the  Mathematical  Statistics  program  under  Task 
SR01 40501  and  Element  61 1 53N  (Job  order  1-1561-108)  and  administered  by  NAVSEA 
Code  05R24. 


INTRODUCTION 


At  the  David  Taylor  Naval  Ship  Research  and  Development  Center  (DTNSRDC), 
random  ship  motion  or  wave  height  data  from  full  scale  trials  or  model  basin 
experiments  is  routinely  subjected  to  various  forms  of  analysis.  Spectral 
analysis  is  performed  to  obtain  autospectra,  cross  spectra  and  response 
amplitude  operators;  and  time  domain  analysis  is  performed  to  obtain  means, 
standard  deviations,  data  point  maxima  and  minima,  data  point  histograms,  peak 
histograms  and  peak-to-trough  histograms.  These  different  types  of  analysis 
reduce  the  time  histories  into  statistical  parameters  that  characterize  the 
test  conditions  and  the  performance  of  the  test  vehicle.  In  particular,  the 
time  history  statistics  that  describe  the  peaks  or  extremes  of  these  random 
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data  are  used  to  predict  the  largest  value  (extreme  value)  that  is  expected  to 
occur  in  a  given  period  of  time.  The  length  of  time  available  for  an  experiment 
is  limited,  so  these  results  must  be  extrapolated  to  longer  time  intervals.  The 
predicted  extremes  in  motion  or  wave  induced  force  are  then  applied  to  the 
design  of  ships  and  other  marine  structures. 

Despite  the  number  of  studies  undertaken  in  this  area,  two  major  problems 
remain.  A  consistent,  easy  to  apply  method  is  required  for  specifying  peak 

occurrences  for  random  processes  with  arbitrary  bandwidths.  And,  second,  a 
general  method  for  accounting  for  statistical  dependence  between  these  events 
is  necessary  when  making  extreme  value  estimates.  This  report  presents  an 
approach  for  solving  these  problems  under  the  assumption  that  the  random 

process  is  GaUssian. 

The  proposed  solution  to  these  problems  is  based  on  applying  a  procedure 
that  produces  the  time  history  of  the  envelope  of  the  actual  test  data;  this 
new  time  history  is  then  analysed  to  produce  the  information  necessary  for 
extreme  value  predictions.  The  envelope  time  history  is  generated  by 
approximating  the  Hilbert  Transform  with  a  digital  filter.  For  a  Gaussian 

random  process,  the  statistics  of  this  envelope  time  history  are  known;  the 
envelope  random  process  is  Rayleigh  distributed  for  wideband  as  well  as 
narrowbanded  time  histories.  The  envelope  time  history  is  then  analysed  to 
obtain  all  standard  statistics  including  a  data  point  histogram  (For 

narrowbanded  time  histories,  this  histogram  is  equivalent  to  the  peak-to-trough 
histogram.)  A  method  developed  for  estimating  or  measuring  the  time  between 
statistically  independent  samples  of  the  envelope  can  then  be  applied.  Once  the 
number  of  statistically  independent  samples  is  known  for  a  given  time,  extreme 
value  predictions  are  then  made. 

For  narrowbanded  random  processes,  each  maximum  and  minimum  does  touch  the 
envelope;  however,  for  random  processes  with  wide  bandwidths,  these  maxima  and 
minima  will  generally  be  lower  than  the  envelope.  In  this  sense,  the  envelope 
approach  can  overpredict  the  extreme  value.  In  the  most  widebanded  simulation, 
the  average  overprediction  was  6.6$.  This  drawback  is  balanced  by  the  method's 
ability  to  determine  if  the  measured  data  is  non-Gaussian  by  performing  a 
Chi-Square  goodness-of-f it  test  on  the  histogram  of  the  envelope  time  history. 
The  envelopes  of  Gaussian  random  processes  with  arbitrary  bandwidths  are  still 
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Rayleigh  distributed.  Since  the  results  of  this  test  are  not  confused  by  the 
bandwidth  of  the  process,  the  characteristics  of  non-Gaussian  processes  are 
more  easily  studied. 

The  method  is  implemented  using  digital  computer  algorithms;  simulated, 
Gaussian  distributed  time  histories  with  varying  bandwidths  are  used  to  develop 
the  procedures.  This  approach  is  verified  by  using  Chi-Square  goodness-of-f it 
tests  to  demonstrate  that  the  envelope  time  history  is  Rayleigh  distributed  and 
that  the  extreme  values  follow  the  distribution  predicted  by  order  statistics. 
This  method  is  then  applied  to  experimentally  obtained  time  histories. 


BACKGROUND 

The  methods  described  in  the  literature  for  making  extreme  value 
predictions  are  based  on  one  or  more  of  the  following  assumptions  about  the 
random  time  history:  The  random  process  is  stationary,  zero  mean  and  Gaussian 
distributed;  its  autospectrum  is  narrowbanded;  and  the  peaks  or  maxima  of  the 
time  history  are  statistically  independent.  The  narrowbanded  assumption  implies 

that  only  one  maximum  and  minimum  (peak-to-trough )  occur  between  two  successive 

1  * 

zero  crossings  of  the  random  process.  The  method  described  by  Longuet-Higgins 
and  by  Cartwright  and  Longuet-Higgins2  is  the  classical  approach.  With  this 
approach,  the  peak-to-trough  amplitudes  are  distributed  according  to  a  Rayleigh 
distribution  for  narrowbanded,  Gaussian  random  processes;  order  statistics 
methods  predict  the  extreme  value  provided  these  amplitudes  are  statistically 
independent  of  each  other.  This  method  requires  only  one  parameter  for  the 
Rayleigh  distribution,  either  the  estimated  mean  square  of  the  peak-to-trough 
amplitudes  or  the  estimated  variance  of  the  time  history.  The  number  of 
peak-to-trough  occurrences  divided  by  the  record  length  determines  the 
encounter  rate;  this  rate  is  then  used  with  order  statistics  methods  to 
estimate  (extrapolate)  the  exposure  time  necessary  for  an  extreme  value  to 
exceed  a  specified  threshold  with  a  given  probability  of  occurrence. 

To  avoid  the  narrowbanded  assumption,  other  approaches  are  described  by 

£ 

A  complete  listing  of  references  is  given  on  page  63. 


3 


Cartwright  and  Longuet-Higgins2  and  by  Ochi^.  These  methods  are  based  on  the 
study  presented  by  Rice14.  Since  a  non-narrowbanded  random  process  can  have  many 
maxima  or  minima  between  successive  zero  crossings,  these  maxima  or  minima  are 
shown  to  be  distributed  according  to  the  distribution  (Rice  distribution)  given 
by  Rice14.  When  applying  order  statistics,  statistical  independence  of  these 
events  is  3till  assumed.  This  method  requires  the  estimation  of  the  first  three 
non-zero  spectral  moments;  these  moments  are  equivalent  to  the  variance  of  the 
random  process  and  the  variance  of  the  first  two  derivatives  of  the  random 
process.  These  moments  are  estimated  from  the  time  history  and  they  are  used  to 
form  a  bandwidth  parameter;  the  Rice  distribution  is  described  by  this 
bandwidth  parameter  and  the  variance  estimate.  The  occurrence  rate  of  these 
events  establishes  the  time  base  for  extreme  value  estimates. 

In  practice,  the  peak-to-trough  approach  proposed  by  Longuet-Higgins1  is 
still  generally  applied  to  experimental  data.  It  is  uncomplicated;  this  method 
only  requires  one  parameter,  the  variance,  to  specify  the  distribution.  The 
average  of  the  highest  one-third  peak-to-trough  amplitudes  is  given  by  four 
times  the  standard  deviation.  Measurement  noise  in  the  test  data  will  add  to 
the  variance  estimate,  and  computation  of  the  peak-to-trough  histogram  mav 

require  a  deadband  zone  such  that  only  the  "significant"  zero  crossing  cycles 
are  used  as  events;  however,  this  approach  is  fairly  tolerant  of  noise.  In 
comparison,  the  method  using  the  Rice  distribution  also  requires  the  variance 
as  well  as  the  bandwidth  parameter.  This  method  is  more  sensitive  to 

measurement  noise  in  the  test  data;  noise  will  change  the  shape  of  the 
histogram  as  well  as  inflate  the  estimated  values  of  the  second  and  fourth 
spectral  moments  that  are  then  used  to  compute  the  bandwidth  parameter. 

Estimation  of  these  spectral  moments  involves  single  and  double  differentiation 
of  the  time  history  (These  calculations  are  usually  performed  in  the  frequency 
domain.);  the  level  of  the  noise  is  increased  by  the  differentiation  process. 

The  peak-to-trough  approach  is  also  applied  to  cases  where  the 

experimental  data  is  non-Gaussian  or  non-narrowbanded.  For  Gaussian  random 
processes  with  moderate  bandwidths,  Tayfun^  proposed  a  peak-to-trough 
probability  density  function  that  is  based  on  sampling  the  Rayleigh  distributed 
envelope.  In  some  cases,  empirical  distributions  such  as  those  given  by  Ochi^ 
or  Forristall^  are  fitted  to  the  measured  peak-to-trough  histograms.  In  other 


cases,  however,  wideband  data  is  simply  assumed  to  be  narrowbanded.  Between  two 
"significant"  zero  crossings,  the  difference  between  the  largest  and  the 
smallest  data  values  is  used  to  form  a  "peak-to-trough"  amplitude.  Attempts  are 
made  to  treat  the  wideband  nature  of  the  data  as  if  it  were  noise.  Deadband 
regions  are  set  to  determine  if  a  zero  crossing  is  "  \ gnif icant" .  Defining  what 
constitutes  a  "significant"  zero  crossing  (or  alternately  a  "significant"  peak 
to  trough  occurrence)  is  a  major  drawback  with  this  approach;  the  required 
criteria  cannot  be  specified  in  a  consistent  manner. 

These  previous  methods  do  not  offer  any  procedures  for  taking  into  account 
the  correlation  between  successive  events  when  making  extreme  value 

Q 

predictions.  Epstein  ,  however,  has  proposed  such  a  procedure  using  a  Markov 

chain  condition.  The  Markov  chain  condition  requires  a  certain  kind  of 

q  1  n 

dependence  between  successive  events.  Recently  Ochi ^  and  Naess  suggest  that 
this  approach  may  account  for  the  effects  of  this  correlation  for  random 
processes  with  moderate  bandwidths;  however,  they  do  not  demonstrate  if  the 
Markov  chain  condition  is  appropriate  for  arbitrary  bandwidth  Gaussian  random 
processes.  Naess10  disagrees  with  Ochi's^  results.  Also,  Epstein®  states  that  a 
clean-cut  analytic  solution  is  expected  for  only  certain  kinds  of  stochastic 
dependence.  This  approach  appears  to  be  subject  to  continuing  development. 

A  procedure  for  estimating  extreme  values  using  the  envelope  associated 
with  the  time  history  is  proposed  by  Naess11;  a  method  for  taking  the  effects 
of  correlation  into  account  is  also  presented.  The  average  rate  at  which  the 
envelope  crosses  the  standard  deviation  level  from  below  is  used  as  part  of  the 
criteria  for  determining  the  effect  of  dependent  samples.  This  procedure  is 
shown  to  be  limited  to  narrowbanded  Gaussian  processes. 

THE  ANALYTIC  ENVELOPE  AND  ITS  DISTRIBUTION 

Regardless  of  the  bandwidth  of  a  Gaussian  distributed  random  process,  the 
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envelope  of  this  process  is  Rayleigh  distributed.  (See  Middleton  ,  pg  351  and 
pg  396;  also,  see  Arens1®  and  Dugundj i 1  tf. )  This  envelope  is  defined  by  using 
the  Hilbert  transform.  In  a  recent  paper  by  Rice1®,  this  envelope  is  called 
the  analytic  envelope. 


The  envelope,  r(t),  is  defined  by 


r(t)  =  [  s(t)2  +  s(t)2  ]1/2 

and,  the  phase  angle,  8,  is  defined  by 

8 ( t )  =  arctan(  s(t)/s(t)  ) 


such  that 


(1) 


(2) 


s(t)  =  r(t)  cos  8(t)  (3) 

where  s(t)  is  the  input  time  history  (either  random  or  deterministic),  and 


s(t)  =  r (t )  sin  8(t)  (4) 

where  s(t)  is  the  Hilbert  Transform  of  s(t).  An  example  of  the  envelope  of  a 
random  time  history  with  moderately  wide  bandwidth  is  given  in  Figure  1.  For  a 
real-valued  function,  s(t),  on  -®  <  t  <  •  ,  the  Hilbert  transform,  s(t),  is 
defined  by 


-  1  ®  s(t) 

s(t)  =  -  /  -  dt  (5) 

TT  -«  t  ~  T 


where  the  principal  value  of  the  integral  is  used.  The  Hilbert  transform 
shifts  all  the  frequency  components  in  the  input  time  history  by  90  degrees.  In 
the  frequency  domain,  the  Hilbert  transform  is  equivalent  to  a  filter  with  the 
following  transfer  function,  H(<d): 


H(w)  =  -  i  sgn(w) 


(6) 


where  i  is  the  square  root  of  minus  one  and  sgn(w)  is  the  signum  function 
tsgn(oj)  =  1  for  u  >  0,  0  for  oi  =  0,  and  -1  for  u  <  0  ]. 

For  a  zero  mean  Gaussian  distributed  random  process,  several  relationships 
exist.  s(t0)  and  s(t0)  are  statistically  independent  at  any  arbitrarily 
selected  time,  t0.  Since  the  envelope  is  the  square  root  of  the  sum  of  two 
statistically  independent  random  variables  at  each  instant  of  time,  the 
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envelope  is  Rayleigh  distributed  with  the  following  probability  density 
function 


p(u) 


-  u2/2 

u  e 


(7) 


and  cumulative  distribution, 


P(u) 


2 

-  u<72 
e 


(8) 


where  u  =  r/o.  o  is  the  standard  deviation  of  s(t).  Also,  s(t)  and  s(t)  have 
the  same  autospectrum  and  autocorrelation  functions.  The  envelope  and  phase 
are  statistically  independent,  and  the  phase  angle  is  uniformly  distributed 
between  +n  and  -it.  The  mean  of  the  envelope  is 

E[  u  ]  =  (  tt/2  )1/2  (9) 


where  E[  .  ]  is  the  expected  value  operator.  The  mean  square  of  the  envelope 
(variance  plus  square  of  the  mean)  is 

E[  u2  ]  =2  (10) 

Since  the  envelope  is  Rayleigh  distributed,  the  average  of  the  highest 

one-third  envelope  samples  is  two  times  the  standard  deviation  of  s(t).  These 
relationships  hold  regardless  of  the  bandwidth  of  the  random  process. 

Except  for  a  narrowband  random  process,  the  peaks  of  the  process  do  not 
necessarily  touch  its  associated  envelope.  See  Figure  1  for  an  example  of  a 
moderately  wide  bandwidth  process.  From  Equation  (1),  the  absolute  value  of  the 
random  process,  |s(t)|,  is  equal  to  the  envelope  whenever  s(t)  is  zero.  Since 
s(t)  and  s(t)  have  the  same  autospectrum  and  hence  the  same  zero  crossing 
statistics,  this  equality  will,  on  the  average,  occur  between  each  zero 
crossing.  From  Equation  (3)t  this  equality  will  occur  whenever  the  phase  is 

either  0  or  -it;  and  since  phase  is  statistically  independent  of  the  envelope,  a 

joint  occurrence  at  an  envelope  maximima  is  therefore  not  expected.  This 

characteristic  of  the  envelope  leads  to  overestimation  of  the  extreme  value  of 


8 
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the  random  process;  the  significance  of  the  overprediction  is  discussed  in  the 
section  titled  Effective  Number  of  Independent  Samples  for  the  Envelopes  from 
the  Simulated  Gaussian  Time  Histories. 

The  Hilbert  transform  is  equivalent  to  a  filter.  Methods  for  approximating 
this  filter  are  described  by  Tayfun^  and  Rabiner  and  Schafer1^.  Given  a  time 
history  that  has  been  digitized  for  analysis  by  a  computer,  an  approximate 
Hilbert  transform  of  this  time  history  can  be  produced.  Thus  a  time  history  of 
the  envelope  can  be  created  and  standard  data  analysis  techniques  (power 
spectra,  correlation,  and  histograms)  can  be  applied  to  this  time  history. 

The  method  used  to  approximate  the  Hilbert  transform  is  given  by  Rabiner 
and  Schafer1 wxth  this  method,  a  non-recursive  digital  filter  (finite  impulse 
response,  FIR)  is  formed,  so  the  Hilbert  transform  of  time  histories  with 
arbitrary  lengths  can  be  calculated  using  a  computer.  From  this  reference,  the 
filter  weights  (impulse  response)  given  in  Table  1  for  a  95  element  filter  were 
used;  these  filter  weights,  h(k),  are  reproduced  in  Table  1.  This  filter  has  a 
peak  approximation  error  of  2.2?  over  the  frequency  interval  from  0.01  to  0.^9 
times  the  sample  rate.  Odd-indexed  filter  weights  are  zero  and  the  last  half  of 
the  filter  wights  are  given  by: 

h(k )  =  -  h(K-l-k)  (11) 

for  k  =  0,  1,  ...  ,  K-1  where  K  is  95.  The  non-recursive  digital  filter  has  the 
following  form: 


-  K-1  K+ 1 

s(JT)  =  l  h(k)  s[(  --y-  *  j  -  k)T]  (12) 
k=0 


where  s(jT)  is  the  filter's  input  and  s(JT)  the  output. 


EXTREME  VALUES  OF  THE  ENVELOPE  PROCESS  AND  ITS  EFFECTIVE  SAMPLE  RATE 

The  extreme  value  of  the  envelope  is  derived  by  applying  order  statistics 
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TABLE  1  - 


HILBERT  TRANSFORM  WEIGHTS  FOR  NON-RECURSIVE  DIGITAL  FILTER 
(FROM  RABINER  AND  SCHAFER16,  TABLE  I,  K=95) 


k 

h(k) 

0 

-.0130099 

2 

-.0045718 

4 

-.0053689 

6 

-.0062800 

8 

-.0072616 

10 

-.0083873 

12 

-.0096455 

14 

-.0110350 

16 

-.0126022 

18 

-.0143770 

20 

-.0163895 

22 

-.0186922 

24 

-.0213465 

26 

-.0244424 

28 

-.0281175 

30 

-.0325665 

32 

-.0380864 

34 

-.0451608 

36 

-.0546233 

38 

-.0680547 

40 

-.0888468 

42 

-.1258168 

44 

-.2112989 

46 

-.6363167 

to  the  sampled  envelope.  Given  that  the  envelope  is  sampled  N  times  such  that 
these  samples  are  statistically  independent,  the  extreme  value  is  defined  as 
the  largest  value  of  the  envelope  that  will  occur  in  these  N  samples.  Once 
these  samples  are  arranged  in  ascending  order,  the  largest  normalized  envelope 
value  in  the  ordered  sample,  u  has  the  following  probability  density 
function: 

6<un)  =  N  p(un)  [  P(un)  ]N_1  (13) 

Since  the  envelope  is  Rayleigh  distributed,  p(u  )  is  given  by  Equation  (7)  and 
P(u  )  is  given  by  Equation  (8).  The  cumulative  distribution  is  given  by: 

u 

G(un)  =  ;n  g(v)  dv  =  {  P(un)  }N  (I'D 

The  expected  or  mean  extreme  value  is  given  by 

00 

%  =  “n  ]  «  fQ  un  «(un}  dun  (15) 

To  this  point,  the  application  of  order  statistics  to  the.  sampled  envelope  is 
equivalent  to  the  treatment  of  order  statistics  to  wave  amplitudes  by 
Longuet-Higgins1 . 

When  the  envelope  of  the  time  history  is  approximated  using  algorithms 
implemented  on  a  digital  computer,  the  sample  rate  is  the  same  for  both  the 
envelope  and  the  time  history.  For  a  time  history  of  fixed  length,  Ng  samples 
are  obtained.  In  general,  these  envelope  samples  will  be  highly  correlated.  A 
measure  of  the  effective  sample  rate  is  required  that  will  produce  Ng 
statistically  independent  samples,  where  Ng  is  the  effective  number  of 
statistically  Independent  samples.  Middleton12  shows  that  uncorrelated  envelope 
samples  are  also  statistically  independent,  so  one  procedure  for  estimating  an 
effective  ‘sample  rate  is  to  compute  the  autocorrelation  function  of  the 
envelope  time  history  and  then  sample  the  envelope  at  time  intervals  such  that 
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this  function  decays  below  a  set  threshold.  The  use  of  this  procedure,  however, 
can  result  in  underestimation  of  extreme  events  because  of  the  reduced  sample 
rate;  sampling  may  not  occur  at  the  peaks  of  the  envelope.  A  approach  that  uses 
all  of  the  available  samples  is  required. 

Two  approaches  are  proposed  for  estimating  or  measuring  the  effective 
sample  rate.  With  either  method,  the  available  time  history  is  first  divided 
into  M  segments  where  each  segment  contains  N  actual  samples.  The  largest 

a 

envelope  value  in  each  segment  is  then  found;  from  one  segment  to  another, 
these  extreme  envelope  values  are  assumed  to  be  statistically  independent.  With 
the  first  approach,  a  histogram  of  these  extreme  envelope  values  is  formed; 
this  histogram  is  then  matched  to  the  hypothetical  distribution  given  by 
Equation  (13)  by  varying  the  number  of  independent  samples,  N.  The  width  of 
each  class  interval  is  selected  to  give  an  equal  number  of  occurrences  for  the 
expected  distribution  (Equation  (14));  the  last  class  interval  has  infinite 
width.  Next,  the  Chi-Square  goodness-of-f it  test  is  applied  to  this  histogram 
using  the  procedures  described  by  Bendat  and  Piersol  The  number  of 
independent  samples  that  minimizes  the  Chi-Square  statistic  is  the  estimate  of 
the  effective  number  of  independent  samples,  Ng .  The  effective  sample  rate  is 
Ne  divided  by  the  segment  length  in  seconds.  This  approach  has  the  advantage 
that  the  measured  distribution  of  extreme  values  is  compared  to  the 
hypothetical  distribution. 

The  second  approach  for  estimating  the  effective  number  of  independent 
samples  is  based  on  matching  the  measured  mean  extreme  value  to  the 
hypothetical  mean  given  by  Equation  (15).  First,  the  mean  is  estimated  for  the 
extreme  envelope  values  from  the  M  segments.  Next,  the  number  of  independent 
samples  that  give  the  same  mean  value  is  found  by  searching  a  look-up  table 
produced  by  numerically  integrating  Equation  (15).  Graphically,  this  procedure 
is  shown  in  Figure  2;  here,  each  of  the  M  segments  contains  N  (actual) 
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samples,  and  the  effective  number  of  independent  samples  is  given  by  Ng.  This 
approach  has  several  advantages.  First,  a  confidence  interval  can  be  placed 
around  the  measured  mean  extreme  value  that  can  then  be  transformed  into  a 
confidence  interval  about  the  effective  number  of  independent  samples.  Second, 
this  procedure  can  be  repeated  after  sampling  the  envelope  at  different  rates 
to  determine  the  effects  of  undersampling  the  envelope  time  history. 
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A  combination  of  these  two  approaches  was  used;  the  effective  number  of 
independent  samples,  Ng>  and  its  confidence  interval  was  found  from  mean  value 
estimates  using  the  second  method,  and  then  the  distribution  of  the  extreme 
values  was  checked  at  this  Ng  value  by  using  the  first  approach.  The 
hypothetical  distribution  was  accepted  at  a  5?  level  of  significance  provided  a 
value  was  found  within  the  confidence  interval  that  also  passed  the 
Chi-Square  goodness-of-f it  test.  90 %  confidence  intervals  were  found  for  these 
mean  value  and  Ng  estimates.  Since  the  mean  of  M  samples  (M  large)  from  any 
distribution  will  approach  a  Gaussian  distribution,  the  procedure  for 
determining  the  mean  value  confidence  interval  given  by  Bendat  and  Piersol^ 
is  used.  For  this  application,  the  true  mean  value  is  known  to  be  located 
within  an  interval  of  width,  w  where: 


s  t. 


M-1 ;0 . 05 


M 


1/2 


(16) 


such  that  is  centered  about  the  estimated  mean  value  where  M  is  the  number 

of  segments,  s  is  the  estimated  standard  deviation,  and  t  is  the  Student  t 

n,o 

distribution.  These  limits  on  the  mean  value  estimate  are  then  used  to  find  the 
confidence  interval  for  the  number  of  statistically  independent  samples.  The 
width  of  this  interval  is  W 


effective  number  of  independent  samples  for  the  envelopes  from 

THE  SIMULATED  GAUSSIAN  TIME  HISTORIES 

These  procedures  for  estimating  the  effective  number  of  independent 
samples  were  applied  to  a  series  of  simulated  Gaussian  time  histories  with 
varying  bandwidths.  From  this  study,  the  ratio  of  the  effective  number  of 
samples  to  the  actual  number  of  samnles  (N  /N  )  is  shown  to  be  related  to  the 
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equivalent  statistical  bandwidth  of  the  sampled  envelope  time  history.  The 
results  from  this  study  can  then  be  applied  to  actual  test  data,  or,  if  the 
test  data  is  sufficiently  long,  these  procedures  for  directly  estimating  this 
ratio  could  be  applied  to  these  test  data. 


The  methods  used  to  produce  these  simulated  time  histories  are  described 
in  Appendix  A.  With  filtered  white  noise  produced  using  a  random  number 
generator,  five  different  time  histories  were  generated  that  had  bandwidths 
varying  from  the  extreme  narrowband  to  wideband.  In  each  case,  the  spectral 
shape  was  a  rectangular  block  centered  at  0.51  hertz;  the  bandwidths  for  each 
time  history  were  0.078,  0.176,  0.293,  0.488  and  0.684  hertz.  These  bandwidths 
correspond  to  runs  2202,  2205,  2203,  2208  and  2206.  The  variance  of  each  time 
history  was  set  to  one  volt.  The  duration  of  these  time  histories  was  10,080 
seconds  at  a  sample  rate  of  five  hertz.  The  envelope  time  histories  for  these 
cases  were  obtained  using  the  Rabiner  and  Schafer1 ^  method  to  produce  the 
Hilbert  transform.  In  Appendix  A,  20  second  long  segments  of  each  time  history 
and  its  envelope  are  reproduced.  Histograms  of  each  time  history  were  obtained 
and  the  Chi-Square  tests  were  run  to  test  the  input  time  history  against  the 
Gaussian  distribution,  and  the  envelope  was  tested  for  the  Rayleigh 
distribution.  In  all  cases,  the  hypothetical  distribution  was  accepted  at  a  5% 
level  of  significance.  Examples  of  the  histograms  are  given.  Also  shown  are 
plots  of  their  autospectra  as  well  as  a  plot  of  each  envelope's  autocorrelation 
function. 

The  procedures  for  finding  the  effective  sample  rates  were  applied  to  each 
of  these  five  envelope  time  histories.  Each  time  history  was  divided  into  60 
intervals  (M=60)  that  were  168  seconds  long.  Six  different  sample  rates  were 
used:  5.0,  2.5,  1.0,  0.5,  0.25  and  0.125  samples  per  second.  Each  segment 

contained  841,  420,  168,  84,  42  or  21  samples  (N  samples)  at  these  respective 
sample  rates.  The  mean  and  standard  deviation  of  the  largest  sample  in  each 
segment  was  estimated;  the  mean  value  estimates  are  given  in  Figure  3.  The 
solid  curve  is  the  expected  mean  numerically  calculated  from  Equation  (15). 
These  results  show  the  affect  of  sampling.  As  the  sample  rate  is  reduced,  the 
estimated  mean  is  almost  constant  (same  peak  values  found)  until  the  curve  of 
the  expected  mean  is  approached.  Here,  the  estimated  mean  reduces  in  value  (the 
affect  of  undersampling),  and  next  the  estimated  mean  follows  the  expected  mean 
as  the  sample  rate  is  reduced  further.  Since  the  mean  is  estimated  from  60 
samples,  a  90$  confidence  interval  is  placed  around  these  mean  estimates.  For 
these  cases,  this  interval  has  a  widtn,  W  ,  from  0.140  to  0.236;  a  nominal 
value  for  each  of  these  cases  is  0.200.  The  width  of  this  interval  is  shown  on 
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Figure  3;  it  describes  the  tolerance  of  these  estimates.  A  summary  of  results 
for  the  five  hertz  sample  rate  is  given  in  Table  2. 

After  equating  these  mean  value  estimates  with  the  look-up  table  obtained 
from  Equation  (15),  the  ratio  of  the  effective  number  of  independent  samples  to 
the  actual  number  of  samples  (N  /n)  was  found.  As  a  function  of  the  estimated 
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equivalent  statistical  bandwidth  of  the  sampled  envelope,  plots  of  these  ratios 
are  given  in  Figure  4.  Regardless  of  the  bandwidth  of  the  original  time 
history,  this  N  /N  ratio  is  shown  to  be  related  only  to  this  equivalent 
statistical  bandwidth  of  the  sampled  envelope. 

The  equivalent  statistical  bandwidth  is  described  in  Appendix  B;  for  our 
application,  this  bandwidth  measure  is  related  to  the  temporal  extent  of  the 
sampled  autocorrelation  function.  The  equivalent  statistical  bandwidth,  B  in 
radians  per  second  is  defined  as: 


w0  2 

[  /  G(io1  du  ] 

0 


/  G(w)  doo 
0 

where  w0  =  tt/T  (one  half  the  sample  rate  in  radians  per  second)  and  G(io)  is  the 
autospectrum  of  the  random  process.  From  Equation  (B.10)  in  Appendix  B,  the 
equivalent  statistical  bandwidth  is  related  to  the  sampled  autocorrelation 
function  of  the  random  process  by 


l  C(kT)2  (18) 

k  =  -“ 

Where  C(kT)  =  R(kT)/R(0)  and  R(kT)  is  the  autocorrelation  function  given  by 
R(kT)  =  E[  x(nT)  x(nT  +  kT)  ]  (19) 


where  x(nT)  is  the  random  process  sampled  at  equal  time  intervals,  T.  The 


TABLE  2  -  SUMMARY  OF  RESULTS  FROM  SIMULATED  TIME  HISTORY  AT 

FIVE  HZ  SAMPLE  RATE 


Run  # 

Mean 

W  /Mean 
m 

Ne 

W  /N„ 
n'  e 

“°/Bst 

VNa 

Chi-Square  Test 

N/Ne 

Test  Stat 

2202 

2.835 

0.082 

34 

0.68 

65.52 

0.040 

. 

4.33 

2205 

3.042 

0.056 

62 

0.52 

35.13 

0.074 

12.33 

2203 

3.291 

0.059 

135 

0.64 

20.80 

0.161 

1.04 

13.67 

2208 

3.450 

0.038 

229 

0.45 

11.76 

0.272 

0.98 

7.00 

2206 

3.565 

0.061 

342 

0.79 

8.40 

0.407 

0.96 

7.33 

TABLE  3  -  AVERAGE  OVERPREDICTION  OF  ENVELOPE  MAXIMA 


Run  No. 

Normalized 

Bandwidth 

Average  Overprediction 
in  Percent 

2202 

0.153 

1.9 

2205 

0.345 

1 . 9 

2203 

0.575 

2.7 

2208 

0.957 

4.3 

2206 

1.341 

6.6 

EFFECTIVE  SAMPLE  RATIO 


sample  period,  T,  is  the  reci  pr  ,u. 
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statistical  bandwidth  was  e.-.tirr;  it  t  y  "  .  •  ■  ,-•••■ 

a  finite  number  of  lags  such  t  h.ii  t-y  ■  ■*rt.  r.  :  ■.  .■ 

short  (decays  to  zero)  compared  to  the  avail  ar-l>-  nimb-v  9  .  t.*  . 

The  equivalent  statist  i>  i;  h.,n  >w:  !'  •.  •  •  .  ; -i. .  * :  '  a  -my 

random  process  as  the  sample  period  is  en.i.ng--  l.  .&.•  •  nr.;.-  ;  .-ric  !,  l,  is 

increased  (reduced  sample  rate),  the  autosjru-t  rum  .f  '  ’.<■  od  random  process 

is  aliased  until  a  limit  is  approached  such  t. hat  h  /U4  1  ;  t  i . i .  imp!  ies  that 
C(kT)  is  one  for  k=0  and  zero  elsewhere,  in  this  case,  uncorrel  >t<-i  successive 
samples  are  used  in  Equation  (19).  For  Gaussian  random  processes  arid  for  the 
Rayleigh  distributed  envelope,  this  also  implies  that  statistically  independent 
samples  are  obtained  at  this  sample  rate.  Next,  as  the  sample  period  is  reduced 
(increased  sample  rate),  the  uutospeetrum  of  the  sampled  run-ion:  process  remains 
unchanged,  so  the  equivalent  statistical  bandwidth  is  constant.  The  sum  of  the 
squared  autocorrelation  samples  (left  hand  side  of  Equation  (18))  is  then 
directly  proportional  to  w0  (one  half  the  sample  rate  in  radians  per  second). 

A  similar  pair  of  asymptotic  limits  are  expected  for  the  Ng/Na  ratios 
given  in  Figure  A.  When  Bst/oj0  =  i,  this  ratio  should  be  one  since 
statistically  independent  samples  are  obtained  for  the  Rayleigh  distributed 
envelope.  For  high  sample  rates,  this  ratio  should  be  inversely  proportional  to 
the  effective  statistical  bandwidth,  since  the  number  of  independent  samples 
becomes  constant. 

The  distribution  of  the  largest  envelope  values  from  each  of  the  r 0 
segments  was  tested  against  the  hypothetical  distribution  given  by  Equation 
(13)  to  check  the  values  estimated  for  N  Using  the  equal  number  of  expected 
occurrences  method,  this  hypothesis  was  tested  using  the  Chi-Square 
goodness-of-f it  test  on  histograms  of  these  extreme  values.  Ten  class  intervals 
were  used,  so  this  test  was  conducted  with  seven  degrees  of  freedom.  As  shown 
in  Table  2,  the  Chi-Square  test  statistic  was  always  less  than  or  equal  to 
14.07,  so,  for  all  five  cases,  this  hypothesis  was  accepted  at  the  S%  level  of 
significance.  This  test  was  performed  using  the  N/Ng  value  given  in  Table  7  for 
the  number  of  statistically  independent  samples.  For  runs  2202,  2203  and  2206, 
the  Chi-Square  test  statistics  plotted  against  N/N^  are  given  in  Figure  6  where 
N  is  allowed  to  vary.  Run  2203  fails  this  test  over  a  wide  range  of  N;  however, 
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Figure  5  -  Chi-Square  Test  Statistic  for  Three  of  the  Simulated 

Gaussian  Time  Histories 
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histograms  with  the  hypothetical  distribution  will  tail  the  Chi -Squ&r< 
five  percent  of  the  time.  For  these  tests  with  variable  N ,  the  i>i!V<'l  opt;  time 
histories  were  sampled  at  five  hertz. 

The  time  history  is  always  bounded  by  its  Rayleigh  distributed  envelope; 
however,  as  seen  in  Figures  A. la  to  A.le  of  Appendix  A,  the  peaks  of  the  *  ime 
history  do  not  always  occur  at  the  peaks  of  its  enveLope.  For  tnis  to  occur, 
the  time  history  maxima  must  coincide  with  an  envelope  maxima  and  a  phase1  angle 
of  zero  or  180  degrees  (See  Equation  (3)-)  Since  the  enveLope  and  phase  angle1 
are  uncorrelated,  such  a  mutual  occurrence  is  unlikely.  This  characteristic  can 
lead  to  overestimation  of  the  extreme  value  reached  by  the  time  history.  The 
size  of  this  overestimation  was  evaluated  using  the  five  simulated  time 
histories.  For  this  study,  each  run  was  again  divided  into  60  segments.  For 
each  segment,  the  largest  envelope  value  and  the  largest  absolute  time  nistu’y 
value  were  found.  The  difference  between  these  two  values  was  divided  by  the 
largest  envelope  value,  and  this  ratio  was  then  averaged  over  the  60  segments 
to  form  a  measure  of  the  average  overprediction.  Expressed  in  percent,  these 
averages  are  given  in  Table  3  where  the  normalized  bandwidth  is  the  bandwidth 
divided  by  the  center  frequency.  For  narrowbanded  cases,  this  average 
overprediction  is  small;  however,  it  increases  with  the  bandwidth  of  the  time 
history.  Even  in  the  most  widebanded  case,  the  predicted  extremes  in  the 
envelope  can  be  used  in  design  applications  for  motion  or  wave  induced  force. 


APPLICATION  TO  EXPERIMENTAL  DATA 

Tnese  procedures  were  applied  to  two  sets  of  experimental  data  obtained 
from  wave  height  measurements.  The  first  data  set  is  a  nyiio  second  segment  from 
the  two  off-shore  platforms  operated  by  the  Naval  Coastal  Systems  Laboratory 
(NCSL)  at  Panama  City,  Florida.  The  second  data  set  consists  of  three  17SF 
second  segments  from  hurricane  Camille  that  occurred  in  1969.  For  each  time 
history.  Appendix  C  contains  plots  of  short  data  segments,  histograms, 
autospectra  and  envelope  autocorrelation  function  plots.  The  histograms  were 
tested  using  the  Chi-Square  goodness-of-f i t  test. 

The  NCSL  wave  height  data  were  taken  at  both  Stages  (off-shore  platforms); 
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Stage  1  is  located  11  miles  (17.7  km)  off-shore  in  100  feet  (30.5  m)  of  water 
and  Stage  2  is  located  2  miles  (3*2  km)  off-shore  in  60  feet  (18. 3  m)  of  water. 
Wave  height  was  measured  using  surface  piercing  wave  gage3.  In  this  report, 
this  data  set  is  labeled  run  2410;  channel  1  is  data  from  Stage  2  and  channel  2 
is  from  Stage  1.  These  data  represent  a  mild  sea  condition;  the  standard 
deviation  was  0.848  feet  at  Stage  1  and  0.526  feet  at  Stage  2.  The  major 
advantage  of  these  data  is  that  the  sea  state  was  reasonably  stationary  over 
the  5940  second  run  length;  tidal  mean  value  shift  was  removed  by  subtracting  a 
best  fit  straight  line  from  these  data. 

The  wave  height  measurements  taken  of  hurricane  Camille  are  documented  by 
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Earle  .  The  envelope  procedures  were  applied  to  three  1752  second  segments 
that  were  taken  at  1015,  1130  and  1500  hours;  these  segments  were  assigned  run 
numbers  2721,  2724  and  2731.  The  Camille  data  is  assumed  to  be  stationary  over 
these  run  lengths.  These  data  represent  a  rough  sea  condition;  the  standard 
deviation  was  5.972  feet  for  run  2721,  5.838  feet  for  run  2724  and  9.966  feet 
for  run  2731 . 

These  wave  data  were  analyzed  using  the  same  procedures  that  were  applied 
to  the  simulated  time  histories.  Each  time  history  was  divided  into  60  segments 
and  the  largest  envelope  value  from  each  segment  was  obtained.  For  the  NCSL 
data,  the  segments  were  99  seconds  long  and  the  following  sample  rates  were 
used:  5.0,  2.5,  1.0,  0.5,  0.25  and  0.125  samples  per  second.  The  number  of 

samples,  Na>  was  495,  247,  99,  49,  24  and  12.  The  Camille  data  was  already 
short,  so  29.2  second  segments  were  used  at  the  following  three  sample  rates: 
2.5,  1.25  and  0.833  samples  per  second.  For  this  data  set,  N  was  73,  36  and 
24.  The  mean  extreme  values  for  these  data  are  given  in  Figure  6;  the  Camille 
data  were  not  undersampled  to  the  point  where  independent  samples  were 
obtained;  however,  both  sets  of  results  show  the  same  trends  as  those  seen  in 
Figure  3>  For  the  NCSL  data,  the  90$  confidence  interval  ranged  from  0.188  to 
0.248,  so  0.220  is  a  nominal  width  for  this  interval.  The  confidence  interval 
for  the  Camille  data  varied  from  0.244  to  0.304;  a  nominal  width  for  this 
interval  is  0.275.  A  summary  of  results  for  the  highest  sample  rates  is  given 
in  Table  4. 

After  equating  mean  extreme  value  estimates,  the  ratio  of  the  effective 
number  of  independent  samples  to  the  actual  number  of  samples  (Ng/Na)  is  given 
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TABLE  4  -  SUMMARY  OF  RESULTS  FROM  MEASURED  WAVE  DATA  AT 

HIGHEST  SAMPLE  RATE 


Run  # 

Mean 

Wm/Mean 

Ne 

w  /nq 

rr  e 

W°/Bst 

VNa 

Chi-Square  Test 

N/N 

e 

Test  Stat. 

2410  Ch  1 

3.376 

0.062 

179 

0.71 

11 . 21 

0.362 

wm 

7.67 

2410  Ch  2 

3.280 

0.076 

130 

0.83 

19.10 

0.263 

KS9 

16.33 

2721 

2.400 

0.126 

11 

0.64 

22.86 

0.151 

Ha 

4.33 

2724 

2.273 

0.135 

8 

0.62 

32.50 

0.110 

0.75 

17.00 

2731 

2.353 

0.129 

10 

0.70 

27.84 

0.137 

0.90 

8.67 

TABLE  5  -  HISTOGRAMS  OF  ENVELOPE  MAXIMA  FOR  SELECTED  MEASURED  WAVE  DATA 


Run  2410 

Ch  2 

Run  2724 

— 

Run  2731 

N 

e 

=  130 

Ne  =  8 

Ne  =  10 

Interval 

#  of 

Interval 

#  of 

Interval 

HQH 

from 

to 

Occur 

from 

to 

Occur 

from 

to 

■ 

0.00 

2.84 

12 

0.00 

1.66 

10 

0.00 

1.78 

10 

2.84 

2.97 

8 

1.66 

1.85 

10 

1.78 

1.95 

9 

2.97 

3.06 

3 

1.85 

1.98 

5 

1.95 

2.09 

4 

3.06 

3.15 

5 

1.98 

2.11 

6 

2.09 

2.21 

7 

3.15 

3.24 

3 

2.11 

2.23 

3 

2.21 

2.33 

4 

3.24 

3.33 

8 

2.23 

2.36 

1 

2.33 

2.45 

4 

3.33 

3.44 

0 

2.36 

2.50 

6 

2.45 

2.59 

4 

3.44 

3.57 

5 

2.50 

2.68 

4 

2.59 

2.76 

4 

3.57 

3.77 

2 

2.68 

2.94 

3 

2.76 

3.02 

4 

3.77 

QD 

14 

2.94 

00 

12 

3.02 

00 

10 

Test 

Stat:  30.00 

Test 

Stat:  19.33 

Test 

Stat :  11.00 

TABLE  6  -  AVERAGE  OVERPREDICTION  OF  ENVELOPE  MAXIMA 
FOR  MEASURED  WAVE  DATA 


Run  # 

Average  Overprediction  in  Percent 

2410  Ch  1 

9.8 

2410  Ch  2 

8.3 

2721 

7.4 

2724 

8.3 

2731 

5.2 

in  Figure  7  as  a  function  of  the  estimated  equivalent  statistical  bandwidth. 
The  solid  curve  is  the  curve  fitted  to  the  simulated  time  history  results  given 
in  Figure  4.  These  results  from  the  wave  data  closely  match  the  results  from 
the  simulated  time  histories.  The  N  /N  ratio  is  again  shown  to  be  related  to 
the  equivalent  statistical  bandwidth  of  the  sampled  envelope. 

The  distribution  of  the  largest  envelope  values  from  each  of  the  60 
segments  was  tested  against  the  hypothetical  distribution  given  by  Equation 
(13).  The  Chi-Square  goodness-of-f i t  test  results  are  given  in  Table  if ;  the 
equal  number  of  expected  occurrences  method  was  used  with  a  total  of  10 
intervals.  For  runs  2410  (channel  2)  and  272*1,  the  Chi-Square  test  statistic  is 
greater  than  1*1.07,  so  the  test  failed  at  the  5%  level  of  significance.  The 
Chi-Square  test  statistic  plotted  against  N/Ng  is  given  in  Figure  8.  In  those 
cases  where  the  test  failed,  larger  numbers  of  occurrences  were  seen  in  the 
first  and  last  intervals;  the  extreme  envelope  values  appear  to  be  scattered 
over  a  wider  range  than  expected  by  the  hypothetical  distribution.  Histograms 
for  these  runs  at  N  =  Ng  are  listed  in  Table  5;  the  histogram  for  run  2731  is 
included  to  give  an  example  of  a  case  that  passed  the  Chi-Square  test  even 
though  it  to  shows  wide  scatter. 

The  average  overprediction  of  the  envelope  maxima  compared  to  the  time 
history  maxima  is  given  in  Table  6  for  each  run.  The  overpredi ct ion  is  several 
percent  larger  than  the  values  seen  for  the  simulated  time  histories. 

When  applying  these  procedures  to  experimental  data,  the  Ng/Na  ratio  for  a 
particular  run  is  the  principal  result  that  is  then  used  to  estimate  extreme 
value  statistics  using  the  methods  of  order  statistics.  Since  N  (actual) 
samples  were  obtained  over  a  known  time  period,  the  effective  number  of 
statistically  independent  samples,  N  ,  over  this  time  period  is  N a  multiplied 
by  this  ratio.  This  time  period  is  known,  so  this  ratio  can  be  converted  to  the 
effective  sample  rate  that  gives  statistically  independent  samples  (Np  divided 
by  this  time  period).  Extrapolation  to  longer  exposure  times  is  accomplished  by 
varying  this  time  interval  and  calculating  a  new  value  for  N  The  expected  or' 
mean  extreme  value  for  this  time  period,  for  example,  is  calculated  using 
Equations  ( 1 3)  and  (15)  where  N  is  equal  to  Ne>  The  Ng/Na  ratio  is  estimated  or 
measured  from  the  available  experimental  data,  so  estimates  of  this  mean 
extreme  value  are  obtained. 
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Figure  8  -  Chi-Square  Test  Statistic  for  Wave  Height  Data 


CONCLUSIONS 


In  the  design  of  ships  and  other  marine  structures,  the  ability  to 
accurately  predict  extremes  in  wave  induced  motions  or  forces  is  limited  by  the 
presence  of  correlated  samples  in  the  measured  data.  This  report  presents  a 
procedure  for  measuring  the  influence  of  these  correlated  samples  in  this 
measured  data;  corrections  are  then  applied  to  these  extreme  value  predictions. 
The  concept  of  an  analytic  envelope  for  arbitrary  bandwidth  Gaussian  random 
processes  is  employed  in  determining  this  correction;  however,  extreme  value 
predictions  are  still  made  within  the  framework  of  standard  order  statistics. 
This  correction  is  shown  to  be  related  to  the  equivalent  statistical  bandwidth 
of  the  random  process's  envelope  time  history.  This  relationship  is 
demonstrated  for  simulated  data  and  for  experimentally  obtained  wave  data. 

The  analytic  envelope  definition  provides  a  consistent,  reasonably  easy  to 
apply  method  for  specifying  peak  occurrences  of  arbitrary  bandwidth  Gaussian 
random  processes.  Since  the  envelope  is  Rayleigh  distributed,  this  approach 
retains  many  of  the  features  found  in  the  peak-to-trough  approach  given  by 
Longuet-Higgins1 .  Digital  filtering  is  required  to  produce  the  envelope  time 
history,  so  the  computer  time  required  to  process  experimental  data  is  a  factor 
in  selecting  this  approach. 

In  a  given  time  history,  the  envelope  maxima  are  generally  larger  than  the 
time  history  maxima.  The  worst  cases  seen  were  an  average  overprediction  of 
6.6?  for  the  simulated  time  histories  and  9.8$  for  the  measured  wave  data. 
Since  the  envelope  is  an  upper  bound  for  a  given  time  history,  the  extreme 
value  predictions  will  be  conservative;  these  predictions  can  be  used  for 
design  applications  for  motion  or  wave  induced  force.  When  extrapolating  to 
longer  exposure  times,  an  important  question  is  the  probability  law  obeyed  by 
the  extreme  values;  the  Chi-Square  test  performed  on  the  largest  envelope 
values  from  each  segment  will  demonstrate  if  the  assumed  distribution  is 
adequate . 


APPENDIX  A 

SIMULATED  GAUSSIAN  TIME  HISTORIES  AND  THEIR  ENVELOPES 


Tne  methods  used  to  produce  the  simulated  time  histories  are  described  in 
this  appendix.  Also  presented  here  are  20  second  long  plots  of  each  time 
history  and  its  envelope.  Histograms  of  each  time  history  were  obtained  and 
Chi-Square  tests  were  run  to  test  the  input  time  history  against  the  Gaussian 
distribution  and  the  envelope  was  tested  for  the  Rayleigh  distribution. 
Examples  of  the  histograms  are  given.  Also  shown  are  plots  of  their  autospectra 
as  well  as  a  plot  of  each  envelope's  autocorrelation  function. 

The  method  for  determining  the  effective  sample  rate  used  in  the  '•’dor 

statistics  estimates  requires  this  series  of  simulated  Gaussian  time  r.ist ~r  i*?n . 

These  time  histories  were  produced  using  algorithms  implemented  on  a  digital 

computer.  Successive  calls  to  the  random  number  generator  produced  bandlimiled 

white  noise;  this  time  history  was  then  filtered  using  a  non-recursi ve  digital 

filter  to  achieve  the  desired  spectral  shape.  The  filter  was  designed  using 
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procedures  similar  to  those  described  by  Helms  .  Since  the  random  number 
generator  produces  uniformly  distributed,  statistically  independent  samples  and 
thirty  of  these  samples  were  averaged  to  produce  one  sample  of  the  band! imi ted 
white  noise,  the  distribution  of  this  noise  and  the  filtered  noise  is 
approximately  Gaussian. 

Five  different  time  histories  were  produced  that  had  bandwidths  varying 
from  the  extreme  narrowband  to  wideband.  In  each  case,  the  spectral  shape  was  a 
rectangular  block  centered  at  0.51  hertz.  The  bandwidths  were  0.073,  0.17b, 
0.293,  0.A88  and  0.681)  hertz.  These  bandwidths  correspond  to  runs  2202,  2205, 
2203,  2208  and  2206.  The  variance  of  each  time  history  was  set  to  one  volt.  The 
duration  of  these  time  histories  was  10,080  seconds  at  a  sample  rate  of  five 
hertz.  The  envelope  time  histories  were  then  obtained  using  the  Rabiner  and 
Schafer^  method  to  produce  the  Hilbert  transform.  Figures  A.  la  through  A.le 
show  a  short  segment  of  each  time  history;  the  envelope  is  plotted  on  each  side 
of  the  axis.  The  segments  shown  in  Figure  A.l  present  the  most  extreme  event 
that  occurred  in  each  time  history.  For  the  narrowbanded  cases,  the  envelope 
tends  to  follow  the  usual  concept  of  a  line  drawn  from  one  time  history  maximum 
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to  the  next;  however,  this  is  not  evident  in  the  wideband  eases.  The  estimated 
or  measured  autospectra  of  these  time  histories  are  given  in  Figure  A. 2.  The 
envelope  autospectra  are  a  maximum  at  or  near  zero  hertz;  they  tend  to  show  the 
triangular  ..hape  that  would  be  expected  from  convolving  the  time  history 
spectrum  with  itself.  The  frequency  range  over  which  they  are  non-zero  appears 
to  be  equal  to  the  range  over  which  the  time  history  spectra  are  non-zero.  The 
estimated  autocorrelation  functions  of  each  envelope  time  history  are  given  in 
Figure  A. 3.  As  expected,  the  narrowbanded  spectra  have  autocorrelations  that  go 
more  slowly  to  zero.  These  autocorrelation  functions  are  also  triangular  in 
shape.  A  summary  of  the  basic  statistics  of  each  time  history  is  given  in  Table 
A .  1 .  In  Table  A.1,  the  normalized  bandwidth  is  the  bandwidth  divided  by  the 
center  frequency  (0.51  hertz).  The  standard  deviation  is  estimated  from  the 
Gaussian  time  history.  The  normalized  maximum  is  the  most  extreme  level 
(absolute  value)  that  occurred  in  that  time  history  divided  by  the  standard 
deviation.  The  normalized  envelope  RMS  is  the  square  root  of  the  variance  plus 
the  mean  squared  for  the  envelope  time  history  divided  by  the  standard 
deviation;  as  seen  from  Equation  (10),  this  ratio  should  be  1.A1A.  In  each 
case,  the  estimated  RMS  is  very  close  to  this  expected  value.  The  normalized 
envelope  maximum  is  the  most  extreme  level  of  the  envelope  time  history  divided 
by  the  standard  deviation.  The  normalized  peak-to-trough  maximum  is  the  most 
extreme  peak-to-trough  of  the  Gaussian  time  history  divided  by  the  standard 
deviation.  The  peak-to-trough  values  are  divided  by  two  to  make  them  comparable 
to  the  envelope  results. 

Data  point  histograms  were  obtained  from  each  time  history  and  its 
envelope.  Using  these  histograms,  the  Chi-Square  goodness-of-f it  test  was 
performed  over  a  range  of  at  least  three  times  the  standard  deviation  of  the 
time  history  to  determine  the  fit  to  the  Gaussian  distribution  for  the  time 
history  and  to  the  Rayleigh  distribution  for  the  envelope.  The  results  of  these 
tests  are  given  in  Table  A. 2;  in  each  case,  the  hypothesis  was  accepted  at  the 
5%  level  of  significance.  The  number  in  parenthesis  under  the  Chi-Square 
statistic  is  the  threshold  for  acceptance  or  rejection  of  the  goodness-of-f it 

test.  The  equal  interval  Chi-Square  goodness-of-f it  procedures  given  by  Bendat 
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and  Piersol  '  were  used.  Since  this  test  requires  statistically  independent 
samples,  the  two  time  histories  were  sampled  at  the  sample  period  given  in 
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Table  A. 2.  The  estimated  autocorrelation  of  each  time  history  has  decayed  well 
into  the  random  noise  level  at  this  time  interval. 

For  runs  2202,  2203  and  2206,  data  point  and  peak-to- trough  histograms  are 
presented  in  Figures  A, 4,  A. 5  and  A. 6;  these  histograms  supplement  the 
Chi-Square  test  results  by  giving  a  visual  fit  to  the  hypothetical 
distributions.  The  data  point  histogram  for  the  time  history  is  plotted  with 
the  Gaussian  distribution,  and  the  data  point  histogram  for  the  envelope  as 
well  as  the  peak-to-trough  histogram  is  plotted  with  the  Rayleigh  distribution. 
To  make  the  scales  comparable  to  the  envelope  data  point  histogram,  the 
peak-to-trough  values  were  divided  by  two  prior  to  computing  their  histogram. 
Every  sample  was  used  to  create  these  histograms;  these  are  not  the  same 
histograms  that  were  used  in  the  Chi-Square  test.  The  Chi-Square  test  was  not 
performed  on  the  peak-to-trough  histograms  since  statistically  independent 
samples  are  required  for  this  test.  The  two  data  point  histograms  consistantly 
match  their  hypothetical  distribution.  As  the  bandwidth  is  increased,  however, 
the  number  of  occurrences  in  the  main  body  of  the  peak-to-trough  histogram  is 
seen  to  increase,  while  the  number  of  occurrences  in  the  tail  begin  to 
decrease;  this  behaviour  is  predicted  by  Tayfun5.  In  these  figures,  the  number 
of  occurrences  in  equal  width  class  intervals  are  plotted  at  each  class 
midpoint . 


TABLE  A. 1 


-  BASIC  STATISTICS  OF  THE  SIMULATED  TIME  HISTORIES  AND 
THEIR  ENVELOPES  (OVER  ENTIRE  RUN) 


Time  History 

;  Envelope 

Peak-to-Trough 
(divided  by  2) 

Run  # 

Band¬ 

width 

Sigma, o 
volts 

Maximum 

RMS  , 

! 

Maximum 

Maximum 

2202 

0.153 

1.015 

5.070 

1.409 

5.080 

5.034 

2205 

0.345 

0.980 

4.199 

1.414 

4.339 

4.007 

2203 

0.575 

0.973 

4.373 

1.411 

4.378 

4.087 

2208 

0.957 

1.016 

3 . 920 

1.412 

4.083 

3.740 

2206 

1 . 341 

0.994 

5.052 

1 . 412 

5.773 

3.692 

TABLE  A. 2  -  CHI-SQUARE  GOODNESS-OF-FIT  TEST  FOR  SIMULATED  TIME 
HISTORIES  AND  THEIR  ENVELOPES 


Run  # 

' 

Sample 

Period 

sec 

Gaussian  Distribution 
(Time  History) 

Rayleigh  Distribution 
(Envelope ) 

Degrees  of 
Freedom 

Chi-Square 
Test  Stat. 

1 

Degrees  of 
Freedom 

| 

Chi-Square 
Test  Stat. 

2202 

16.0 

18 

21,47  * 

8 

8.83 

(28.87) 

(15.51) 

2205 

13.0 

17 

20.73 

10 

13.61 

(27.59) 

(18.31) 

2203 

10.0 

18 

15.18 

8 

7.24 

(28.87) 

(15.51) 

2208 

6.0 

19 

18.27 

8 

11.49 

(30.14) 

(15.51) 

2206 

4.0 

29 

27.55 

9 

5.80 

(42.56) 

(16.92) 

* 

Rejection  threshold 

for  Chi-Square  Test 

34 
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Figure  A. 2. a  -  Autospectra  of  Simulated  Gaussian  Time  Histories 
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APPENDIX  B 

RELATIONSHIP  BETWEEN  EQUIVALENT  STATISTICAL  BANDWIDTH  AND 
THE  AUTOCORRELATION  FUNCTION 


The  definition  of  the  equivalent  statistical  bandwidth  of  a  filter  is 
given  by  Bendat  and  Piersol1"^;  their  definition  is  now  extended  to  the 
equivalent  statistical  bandwidth  of  a  random  process  with  any  distribution. 
From  Bendat  and  Piersol,  this  "is  the  bandwidth  of  a  hypothetical  rectangular 
filter  which  would  pass  a  signal  with  the  same  mean  square  value  statistical 
error  as  the  actual  filter  when  the  input  is  white  noise".  The  mean  square 
value  statistical  error  relationships  derived  by  Bendat  and  Piersol  assume  a 
Gaussian  distributed  random  process;  however,  this  bandwidth  definition  has 
meaning  for  random  processes  other  than  Gaussian.  The  equivalent  statistical 
bandwidth,  B  in  radians  per  second  of  a  random  process  is  defined  as: 

a)0  2 

[  /  G(w)  dw  ] 

0 


/  G(u>)  dw 

0 

where  w0  =  ir/T  (one  half  the  sample  rate  in  radians  per  second)  and  G(uj)  is  the 
autospectrum  of  the  random  process.  Since  this  random  process  is  sampled,  it 
is  bandlimited  (G(w)  =0);  it  has  nonzero  values  only  in  the  frequency  interval 
-w0  ^  u  i  u,,  The  relationship  between  the  equivalent  statistical  bandwidth  and 
the  autocorrelation  function  of  the  random  process  is  now  derived;  this 
derivation  is  independent  of  the  distribution  of  the  random  process.  It  is 
based  on  sampled  random  processes;  however,  this  derivation  can  be  extended  to 
continuous  processes. 

1  Q 

Balakri shnan  7  proved  the  following  relationships  between  autocorrelation 
and  the  autospectrum  for  sampled,  stationary  random  processes.  The 
autocorrelation  function,  R(kT),  is  given  by 

=  E[  x(nT)  x(nT  +  |<T)  ] 

H3 


R(kT) 


( B.  2) 


where  x(nT)  is  the  sampled  time  history.  The  Fourier  transform  pair  are 


R(kT) 


1  <d0  ikTw 

/  G (id )  e  dio 

2  -<d0 


and 


(B.3) 


G(w) 


-ikTw 

R(kT)  e 


(B.  4) 


2 

The  top  integral  in  Equation  (B.1)  is  the  variance  of  the  random  process,  o  , 
the  zeroth  autocorrelation  lag: 


1  “o 

»  R(0)  =  /  G(w)  dm  /  G(id)  dm  (B.5) 

2  -io0  0 


The  lower 

integral 

in  Equation 

(B.1) 

is  found 

by  substituting 

Equation 

for  G (u>) : 

<d0 

1 

QO  00 

(Do 

-iu)T(k+l ) 

/  G(<d)2 

dw  = - 

l  l 

RUT) 

RUT) 

; 

e  dm 

(B.  6) 
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M 

o 
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CM 
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> 

"(Do 
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(D0 

-iuT(k+l) 

2u0  , 

k  +  1 

=  0 

/  e 

da) 

=  1 

( B.  7) 

-<d0 

0, 

otherwi se 

then 

(a)q 

/ 

G(id)2  dm 

1 

00 

l 

R( kT) 2 

(B.8) 

0 

(D0 

k  =  - 

00 

Substituting  Equations  (B.5)  and  (B.8)  into  Equation  (B.1)  gives  the  relation¬ 
ship  between  the  equivalent  statistical  bandwidth  and  the  autocorrelation 
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function  of  the  random  process: 


or 


st 


R(0)2 

w0  - 

00 

l  R(kT)2 

k  =  -“ 


“o 


B 


st 


I  C (kT) 2 
k=~® 


Where  C(kT)  is  the  normalized  autocorrelat ion  function, 


(B.  9) 


(B. 10) 


R  (kT ) 

C  (kT)  = -  (B.  1 1) 

R(  0) 

Equation  (B.9)  gives  the  relationship  between  the  equivalent  statistical 
bandwidth  and  the  autocorrelation  function.  At  a  fixed  sample  rate,  the 
equivalent  statistical  bandwidth  is  inversely  proportional  to  the  sum  of  the 
squared  autocorrelation  samples.  In  general,  the  statistical  bandwidth  is  a 
measure  of  the  extent  (in  time)  of  the  autocorrelation  function;  it  summarizes 
this  feature  of  the  autocorrelation. 


APPENDIX  C 

MEASURED  WAVE  HE I  CUT  DATA 


The  experimental  data  were  obtained  from  two  sets  of  wave  height 
measurements.  The  first  data  set  is  a  5940  second  segment  from  the  two 
off-shore  platforms  operated  by  the  Naval  Coastal  Systems  Laboratory  (NCSL)  at 
Panama  City,  Florida.  The  second  data  set  consists  of  three  1752  second 
segments  from  hurricane  Camille  that  occurred  in  1969.  This  appendix  contains 
plots  of  short  data  segments,  histograms,  autospectra  and  envelope 
autocorrelation  function  plots  for  each  data  run.  Chi-Square  goodness-of-f it 
tests  were  performed  on  the  histograms. 

The  NCSL  wave  height  data  were  taken  at  both  Stages  (off-shore  platforms); 
Stage  1  is  located  11  miles  (17.7  km)  off-shore  in  100  feet  (30.5  m)  of  water 
and  Stage  2  is  located  2  miles  (3.2  km)  off-shore  in  60  feet  ( 1 8. 3  m)  of  water. 
Wave  height  was  measured  using  surface  piercing  wave  gages.  In  this  report, 
this  data  set  is  labeled  run  2410;  channel  1  is  data  from  Stage  2  and  channel  2 
is  from  Stage  1.  These  two  wave  height  records  are  reasonably  stationary  over 
the  5940  second  run  length;  tidal  mean  value  shift  was  removed  by  subtracting  a 
best  fit  straight  line  from  these  data.  The  hurricane  Camille  wave  height  data 
are  documented  by  Earle20.  The  three  1752  second  long  data  segments  selected 
for  anlaysis  were  taken  at  1015,  1130  and  1500  hours;  these  segments  were 
assigned  run  numbers  2721,  2724  and  2731.  The  Camille  data  is  assumed  to  be 
stationary  over  these  run  lengths. 

Both  sets  of  wave  data  were  processed  to  obtain  the  envelope  time  history, 
autospectra,  envelope  autocorrelation  functions,  data  point  histograms  and 
peak-to-trough  histograms,  as  well  as  basic  statistics.  The  envelope  time 
histories  were  then  obtained  using  the  Rabiner  and  Schafer ^  method  to  produce 
the  Hilbert  transform.  A  summary  of  the  basic  statistics  of  each  time  history 
is  given  in  Table  C.l.  The  standard  deviation  is  estimated  from  the  wave  time 
history.  The  normalized  maxima  is  the  most  extreme  level  (absolute  value)  that 
occurred  in  that  time  history  divided  by  the  standard  deviation.  The  normalized 
envelope  RMS  is  the  square  root  of  the  variance  plus  the  mean  squared  for  the 
envelope  time  history  divided  by  the  standard  deviation;  as  seen  from  Equation 
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(10),  this  ratio  should  be  1.41 4.  In  each  case,  the  estimated  RMS  is  very  close 
to  this  expected  value.  The  normalized  envelope  maxima  is  the  most  extreme 
level  of  the  envelope  time  history  divided  by  the  standard  deviation.  The 
normalized  peak-to- trough  iir,.-.ima  is  the  most  extreme  peak-to-trough  of  the  wave 
time  history  divided  by  the  standard  deviation. 

Figures  C.1  through  C.12  present  time  history  segments,  autospectra, 
histograms  and  autocorrelation  function  estimates  from  these  time  histories. 
These  figures  are  organized  by  wave  data  set  and  run  number  in  the  following 
manner:  For  run  2410,  20  second  segments  of  each  time  history  are  given  in 
Figure  C.1  for  channel  1  and  Figure  C.4  for  channel  2;  for  the  Camille  wave 
aata,  60  second  segments  are  given  in  Figures  C.8  for  run  2721,  Figure  C.11  for 
run  2724  and  Figure  C.1 4  for  run  2731.  Each  segment  brackets  the  most  extreme 
event  that  occurred  during  that  run;  the  envelope  is  plotted  on  each  side  of 
the  axis.  For  the  corresponding  runs,  wave  and  envelope  autospectrum  estimates 
are  plotted  in  Figures  C.2,  C.5,  C.9,  C.12  and  C.15.  Similarly,  the  data  point 
histograms  and  peak-to-trough  histograms  are  plotted  in  Figures  C.3,  C.6,  C.10, 
C.1 3  and  C.1 6.  Autocorrelation  estimates  for  the  envelope  time  history  are 

given  in  Figure  C.7  f or  channels  1  and  2  of  run  2410;  the  autocorrelation 

estimates  for  runs  2721,  2724  and  2731  are  given  in  Figure  C.1 7. 

The  envelope  autospectra  given  in  Figures  C.2,  C.5,  C.9,  C.12  and  C.15 
show  the  general  features  of  the  envelope  autospectra  for  the  simulated  data 
presented  in  Appendix  A.  Similarly,  the  autocorrelation  functions  given  in 

Figures  C.7  and  C.1 7  share  the  general  features  of  the  simulated  data 
autocorrelation  functions.  The  wave  autospectrum  for  run  2410,  Channel  1  is 
bimodal;  however,  this  feature  does  not  seem  to  affect  its  envelope  auto 

spectrum  or  its  autocorrelation  function  estimate. 

Data  point  and  peak-to-trough  histograms  are  presented  in  Figures  C.3, 
C.6,  C.10,  C .  1 3  and  C . 1 6 ;  these  histograms  supplement  the  Chi-Square  test 
results  presented  below  by  giving  a  visual  fit  to  the  hypothetical 
distributions,  me  data  point  histograms  for  the  time  histories  (Figures  C.3a, 
C.6a,  C.lOa,  C.13a  and  C.l6a)  are  plotted  with  the  Gaussian  distribution,  and 
the  data  point  histograms  for  the  envelope  (Figures  C.3b,  C.6b,  C.lOb,  C .  1 3b 
and  C.l6b)  as  well  as  the  peak-to-trough  histograms  (Figures  C.3c,  C.6c,  C.IOc, 
C . 1 3c  and  C.l6c)  are  plotted  with  the  Rayleigh  distribution.  Every  sample  was 
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used  to  create  these  histograms;  these  are  not  the  same  histograms  that  were 
used  in  the  Chi-Square  tests  described  below.  Except  for  the  peak-to-trough 
histograms,  a  good  qualitative  match  to  their  hypothetical  distribution  is 
seen.  The  peak-to-trough  histograms  appear  to  have  more  occurrences  near  the 
peak  of  the  hypothetical  distribution,  and  less  occurrences  in  the  tail;  this 
feature  was  more  obvious  with  the  widebanded  simulated  time  histories.  In  every 
case,  the  wave  data  point  histograms  have  occurrences  at  a  higher  level  on  the 
positive  side  of  the  distribution  than  they  do  on  the  negative  side;  non-linear 
wave  effect  probably  account  for  this  feature;  howeve.  ,  except  for  channel  1  of 
run  2410,  the  other  wave  data  passed  the  Chi-Square  test  for  the  Gaussian 

distribution. 

An  additional  set  of  data  point  histograms  was  obtained  from  each  time 
history  and  its  envelope.  Using  these  histograms,  the  Chi-Square 

goodness-of-f it  test  was  performed  over  a  range  of  at  least  three  times  the 
standard  deviation  of  the  time  history  for  channels  1  and  2  of  run  2410  and 
over  a  range  of  at  least  two  times  the  standard  deviation  for  runs  2721,  2724 
and  2731  .  These  tests  were  performed  to  determine  the  fit  to  the  Gaussian 
distribution  for  the  wave  data  and  to  the  Rayleigh  distribution  for  the 

envelope.  The  Chi-Square  test  was  not  performed  on  the  peak-to-trough 

histograms  since  statistically  independent  samples  are  required  for  this  test. 
The  results  of  these  tests  are  given  in  Table  C.2;  except  for  run  2410  channel 
1,  the  hypothesis  was  accepted  at  the  5J  level  of  significance.  The  number  in 
parenthesis  under  the  Chi-Square  statistic  is  the  threshold  for  acceptance  or 
rejection  of  the  goodness-of-f it  test.  The  equal  interval  Chi-Square 
goodness-of-f it  procedures  given  by  Bendat  and  Piersol  '  were  used.  Since  this 
test  requires  statistically  independent  samples,  the  two  time  histories  were 
sampled  at  the  sample  period  given  in  Table  C.2.  The  estimated  autocorrelation 
of  each  time  history  has  decayed  well  into  the  random  noise  level  at  this  time 
interval . 
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TABLE  C . 1 


-  BASIC  STATISTICS  OF  THE  MEASURED  WAVE  DATA  AND  THEIR 
ENVELOPES  (OVER  ENTIRE  RUN) 


Run  # 

Time  History 

Envelope 

Peak-to-Trough 
(divided  by  2) 

Sigma, a 
feet 

Maximum 

RMS 

Maximum 

Maximum 

2410  Ch  1 

0.526 

4.643 

1.411 

4.738 

3.379 

2410  Ch  2 

0.848 

4.476 

1.417 

4.509 

3.445 

2721 

5.972 

3.674 

1.418 

3.878 

2.877 

2724 

5.838 

4.722 

1.419 

4.764 

4.374 

2731 

9.966 

4.827 

1.418 

4.849 

3.692 

TABLE  C. 2  -  CHI-SQUARE  GOODNESS-OF-FIT  TEST  FOR  MEASURED  WAVE  DATA 

AND  THEIR  ENVELOPES 


Run  # 

Sample 

Period 

sec 

2410  Ch  1 

15.0 

2410  Ch2 

15.0 

2721 

30.0 

2724 

30.0 

2731 

30.0 

Gaussian  Distribution 
(Time  History) 


Degrees  of  Chi-Square 
Freedom  Test  Stat. 


30.34 


(27.59) 

11.89 

(25.00) 

4.32 

(12.59) 

8.53 

(14.07) 

8.12 

(14.07) 


Rayleigh  Distribution 
(Envelope) 


Degrees  of  Chi-Square 
Freedom  Test  Stat. 


15.74 

(16.92) 

3.06 

(14.07) 

7.58 

(9.49) 

5.40 

(11.07) 

1.14 

(9.49) 


Rejection  threshold  for  Chi-Square  Test 
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Figure  C.13  -  Histograms  for  Run  2724 


AUTOSPECTRUM  (feet2  sec) 


Figure  C.14  -  Time  Hist 


1800.0 


900.0 


0.0 


0.0 


0.1 


FREQ1 

Figure  C.15  -  Autospectra  of 

for  Run 


300  + 


o  Measured 
x  Gaussian 


NORMALIZED  ENVELOPE  AMPLITUDE 


/ 


•  : 


o 

X 


\ 


a 


V  *• 


‘pS  i 


i  i  » "  i  i  i  "i 


o.o 


Measured 

Rayleigh 


NORMALIZED  PEAK-TO-TROUGH  AMPLITUDE 


Figure  C.16  -  Histograms  for  Run  2731 
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Figure  C.17  -  Autocorrelation  Functions  for  Envelopes  of 
Runs  2721.  2724  and  2731 


300  + 


o  Measured 
x  Gaussian 


Figure  C.16  -  Histograms  for  Run  2731 


REFERENCES 


1.  Longuet-Higgins ,  M.S.,  "On  the  Statistical  Distribution  of  the  Heights 
of  Sea  Waves,"  Journal  of  Marine  Research,  Vol .  11,  No.  3,  p.  245  (1952). 

2.  Cartwright,  D.E.  and  Longuet-Higgins,  M.S.,  "The  Statistical 
Distribution  of  the  Maxima  of  a  Random  Process,"  Proc.  Roy.  Soc.  London,  Ser. 

A,  237,  p.  21 2  (1956) . 

3.  Ochi ,  M.K.,  "On  Prediction  of  Extreme  Values,"  Journal  of  Ship 
Research,  Vol.  17,  No.  1 ,  p.  29  (Mar  1973) - 

4.  Rice,  S.O.,  "Mathematical  Analysis  of  Random  Noise,"  Bell  System 
Technical  Journal,  Vol.  23  0944)  and  Vol.  24  0945). 

5.  Tayfun,  M.A.,  "Distribution  of  Crest-to-Trough  Wave  Heights,"  J.  Wat. 
Ways  Port  Coastal  Ocean  Div.  Am.  Soc.  civ.  Engrs.,  107(WW3),  p.  149  (Aug  1981). 

6.  Ochi,  M.K.,  "Generalization  of  Rayleigh  Probability  Distribution  and 
Its  Application,"  Journal  of  Ship  Research,  Vol.  22,  No.  4,  p.  259  (1978). 

7.  Forristall,  G.Z.,  "On  the  Statistical  Distribution  of  Wave  Heights  in 
a  Storm,"  Journal  of  Geophysical  Research,  Vol.  83,  No.  05,  p.  2353  (1978). 

8.  Epstein,  B. ,  "The  Distribution  of  Extreme  Values  in  Sample  whose 
Members  are  Subject  to  a  Markov  Chain  Condition,"  Annals  of  Mathematical 
Statistics,  Vol.  20,  p.  590  (1949). 

9.  Ochi,  M.K.,  "Extreme  Values  of  Waves  and  Ship  Responses  Subject  to  the 
Markov  Chain  Condition,"  Journal  of  Ship  Research,  Vol.  23i  No.  3,  p.  188  (Sept 
1979) . 


10.  Naess,  A.,  "Extreme  Values  of  a  Stochastic  Process  Whose  Peak  Values 
are  Subject  to  the  Markov  Chain  Condition,"  Norwegian  Maritime  Research,  No.  1, 
Vol.  11,  p.  16  (  1983)  . 

11.  Naess,  A.,  "Extreme  Value  Estimates  Based  on  the  Envelope  Concept," 
Applied  Ocean  Research,  Vol.  4,  No.  3,  p.  181  (1982). 

12.  Middleton,  D.,  An  Introduction  to  Statistical  Communication  Theory, 
McGraw  Hill,  N.Y.,  N.Y.  (I960) . 

13.  Arens,  R.,  "Complex  Processes  for  Envelopes  of  Normal  Noise,"  IRE 
Transactions  on  Information  Theory,  Vol.  IT-3,  p.  204  (Sept  1957). 

14.  Dugundji,  J.,  "Envelopes  and  Pre-Envelopes  of  Real  Waveforms,"  IRE 
Transactions  on  Information  Theory,  Vol.  IT-4,  p.  53  (Mar  1958). 

15.  Rice,  S.O.,  "Envelopes  of  Narrow-Banded  Signals,"  Proceedings  of  the 
IEEE,  Vol.  70,  No.  7,  p.  692  (July  1982). 


63 


A 


REFERENCES  (Cont.) 


16.  Rabiner,  L.R .  and  Schafer,  R.W.,  "On  the  Behavior  of  Minimax  FIR 
Digital  Hilbert  Transformers,"  The  Bell  System  Technical  Journal,  Vol .  53.  No. 
2,  p.  363  (Feb  1974). 

17.  Bendat,  J.  and  A.  Piersol,  Random  Data;  Analysis  and  Measurement 
Procedures,  Wiley-Interscience,  N.Y.,  N.Y.  (1971). 

18.  Helms,  H.D.,  "Nonrecursive  Digital  Filters:  Design  Methods  for 
achieving  Specifications  on  Frequency  Response,"  IEEE  Trans.  Audio  and 
Electroacoustics,  Vol.  AU-16,  No.  3,  p.  336  (Sept  1969). 

19.  Balakri shnan ,  A.V.,  "A  Note  on  the  Sampling  Principle  for  Continuous 
Signals,"  IRE  Trans.  Information  Theory,  vol.  IT-3,  p.  143  (June  1957). 

20.  Earle,  M.D.,  "Extreme  Wave  Conditions  During  Hurricane  Camille," 
Journal  of  Geophysical  Research,  Vol.  80,  No.  3,  p.377  (Jan  1975). 


64 


INITIAL  DISTRIBUTION 


Copies  CENTER  DISTRIBUTION 


6 

NAVSEA 

Copies 

Code 

Name 

5  SEA  05R24/Sejd 

1  SEA  55W3/Comstock 

1 

1600 

H.  R.  Chaplin 

12 

DTIC 

1 

1604 

F.  W.  Wilson 

1 

U.  of  Florida/M.  K.  Ochi 

1 

1630 

A.  G.  Ford 

1 

City  College/W.  J.  Pierson 

1 

1630 

W.  M.  Richardson 

1 

1700 

W.  W.  Murray 

CENTER 

DISTRIBUTION 

1 

1730 

A.  B.  Stavovy 

Copies 

Code 

Name 

1 

1730.6 

W.  H.  Buckley 

1 

012.3 

D .  D .  Moran 

1 

2700 

R.  C.  Allen 

1 

1100 

W.  C.  Dietz 

1 

2800 

J.  R.  Belt 

1 

1173 

Z.  G.  Wachnik 

1 

2900 

G.  G.  Switzer 

1 

1500 

W.  B.  Morgan 

1 

2950 

R.  T.  Schwartz 

1 

1504 

V.  J.  Monacella 

1 

2960 

D.  J.  Plumpe 

1 

1504 

J.  A.  Fein 

10 

2960 

R.  D.  Pierce 

1 

1513 

W.  E.  Cummins 

1 

522.1 

Unclassified  Lib 

1 

1522 

M.  B.  Wilson 

1 

522.2 

Unclassified  Lib 

1 

1560 

D.  3.  Cieslowski 

1 

1561 

G.  G.  Cox 

4 

1561 

J.  F.  0'Dea 

1 

1561 

S.  L.  Bales 

4 

1562 

M.  J.  Davis 

1 

1563 

J.  E.  Russ 

DTNSRDC  ISSUES  THREE  TYPES  OF  REPORTS 

1.  DTNSRDC  REPORTS.  A  FORMAL  SERIES,  CONTAIN  INFORMATION  OF  PERMANENT  TECH 
NICAL  VALUE.  THEY  CARRY  A  CONSECUTIVE  NUMERICAL  IDENTIFICATION  REGARDLESS  OF 
THEIR  CLASSIFICATION  OR  THE  ORIGINATING  DEPARTMENT. 

2.  DEPARTMENTAL  REPORTS.  A  SEMIFORMAL  SERIES,  CONTAIN  INFORMATION  OF  A  PRELIM 
INARY.  TEMPORARY,  OR  PROPRIETARY  NATURE  OR  OF  LIMITED  INTEREST  OR  SIGNIFICANCE. 
THEY  CARRY  A  DEPARTMENTAL  ALPHANUMERICAL  IDENTIFICATION. 

3.  TECHNICAL  MEMORANDA.  AN  INFORMAL  SERIES.  CONTAIN  TECHNICAL  DOCUMENTATION 
OF  LIMITED  USE  AND  INTEREST.  THEY  ARE  PRIMARILY  WORKING  PAPERS  INTENDED  FOR  IN 
TERNAL  USE.  THEY  CARRY  AN  IDENTIFYING  NUMBER  WHICH  INDICATES  THEIR  TYPE  AND  THE 
NUMERICAL  CODE  OF  THE  ORIGINATING  DEPARTMENT.  ANY  DISTRIBUTION  OUTSIDE  DTNSRDC 
MUST  BE  APPROVED  BY  THE  HEAD  OF  THE  ORIGINATING  DEPARTMENT  ON  A  CASE-BY-CASE 
BASIS. 


